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MODELING  OF  CURRENT  FEATURES  IN  GULF  STREAM  SAR  IMAGERY 


It  IntrQdwCtPPW 

Interpretation  of  ocean  radar  imagery  from  regions  involving  strong  current  flows,  such 
as  the  Gulf  Stream,  is  possible  only  if  the  link  between  observed  radar  features  and  the 
underlying  depth  and  surface  current  structure  is  sufficiently  understood.  Because  the  radar 
effectively  maps  the  surface  waves,  this  link  requires  knowledge  of  the  effects  of  wave- 
current/depth  interactions  on  the  wind-wave  field.  These  effects  must  be  understood  for  a  wide 
range  of  wavelengths;  waves  on  the  1  m  or  larger  scales  lead  to  tilt  modulation  effects  in  the 
radar  image  and  cm-scale  waves  affect  the  radar  return  as  Bragg  scatterers.  These  modulations 
in  the  wave  field  often  lead  to  sharp  features  in  radar  imagery  even  for  relatively  mild 
current/depth  variations.  These  features  are  produced  by  both  advective  and  refractive  effects 
which  can  alter  the  size  and  propagation  direction  of  the  waves.  An  obvious  example  of 
refractive  effects  on  waves  is  evident  in  coastal  regions  where  a  sloping  beach  tends  to  orient 
the  propagation  direction  of  waves  toward  the  shore.  Waves  in  coastal  areas  are  also 
significantly  affected  by  current-depth  variations  in  tidal  inlets  or  near  the  mouths  of  rivers.  The 
refractive  effect  of  the  gulf  stream  cunent  on  the  propagation  of  long  transatlantic  ocean  waves 
is  also  easily  observed  and  well  understood.  Much  less  is  known  about  the  wave-current 
interaction  for  shorter  waves  (1  cm  to  100  m),  which  often  is  a  dominant  factor  in  the  formation 
of  high  intensity  variations  in  the  radar  return  that  are  frequently  observed  in  real  and  synthetic 
aperture  radar  imagery. 

To  investigate  current  depth  features  in  the  ocean,  NRL  has  developed  a  Time-dependent 
Ocean  Wave  (NTOW)  model  which  allows  for  general  time  varying  depth,  current,  and  wind 
fields.  The  model  uses  conservation  of  wave  action  to  compute  the  propagation  of  a  statistical 
wind  wave  field.  The  action  density  formalism  is  presented  in  section  2,  and  the  NTOW  model 
impiementation  of  the  solution  of  the  action  density  conservation  equation  using  a  finite  step 
scheme  is  described  in  section  3.  Tests  of  the  NTOW  model  are  presented  in  section  4.  In 
section  5  we  present  NTOW  predictions  for  the  wind  wave  field  and  radar  return  for  a  current 
model  of  the  "Rip"  feature  which  appeared  in  NRL  high  resolution  Rcol  Aperture  Radar  (RAR) 
imagery  of  a  boundary  region  of  the  Gulf  Stream  off  Cape  Hatteras.  Section  6  summarizes. 


Manuscript  approved  September  17,  1993. 
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2.1.  Dispersion  Relation 


For  surface  gravity-capillary  waves  in  the  linear  approximation  in  the  absence  of  currents 
the  dispersion  relation  relating  the  wavenumber  k  and  the  intrinsic  frequency  Uq  is 

»  gk  tanhJbf  \  *  —  ,  (1) 

Pg 

where  g  is  the  acceleration  of  gravity,  d  is  the  water  dq)th,  T  is  the  surface  tension,  and  p  is  the 
water  density.  The  wave  crests  propagate  at  the  phase  speed  c  given  by  the  ratio. 


whereas  the  wave  energy  propagates  at  the  group  speed  Cg  given  by 

C  .  .  fl  *  .  1  1  . 

^  dk  2  sinh2itd  i  +  ^ 


When  non-zero  currents  are  present  the  frequency  of  the  waves  observed  in  the  rest  frame 
will  be  Doppler  shifted.  The  frequency  in  the  rest  frame  is  related  to  the  local  intrinsic 
frequency  (in  the  moving  frame)  by, 

(j)  =  Ug  +  k  •  U.  (4) 

It  is  interesting  to  observe  that  for  sufficiently  large  currents  the  second  term  in  Eq.  4  can 
become  comparable  to  Uo-  In  this  case,  the  frequency  of  the  waves  opposing  the  current  and 
their  energy  propagation  speed,  Cg+fc'U,  can  vani^  leading  to  the  phenomenon  of  wave 
blocking. 


2.2.  Wave  Action  Density  Equation 


For  slowly  varying  mean  fields  the  evolution  of  the  wave  spectrum  is  most  easily 
modeled  in  terms  of  the  wave  action  density  formalism  [cf.  G.B.  Whitham,  1967;  F.P. 
Bretherton  and  C.J.R.  Garrett,  1968].  The  wave  field  is  regarded  as  a  superposition  of  waves, 
such  that  the  free  surface  displacement  f  (x,/)  can  be  represented  by  a  Fourier-Stiel^es  sum, 

e(x.o  =  (s) 

where  ij*  is  a  random  Fourier  amplipide.  For  a  homogeneous,  stationary  wave  field  rjk  satisfies 
<11,1,;,  >  -  •  F(k)  ^  ,  (6) 

where  Ak  is  the  two  dimensional  wavenumber  increment,  E(k)  is  the  energy  density  spectrum, 
and  F(k)  is  the  variance  density  spectrum  satisfying 

=  <|2>.  (7) 


In  the  absence  of  depth  or  current  variations  the  variance  spectrum  is  conserved  along 
ray  paths.  However,  when  current  or  depth  variations  are  present  F  is  no  longer  conserved  due 
to  the  variations  in  the  frequency  along  these  ray  paths.  The  wave  action  density  defined  by 


N(k-^,t) 


E(k;x,0  ^ 

k 


(8) 


however,  remains  constant  along  ray  paths  even  for  variable  depths  and  currents  when  no  energy 
sources  or  sinks  are  present.  The  wave  action  density  conservation  equation  for  a  source  S  can 
be  written  as 
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(9) 


M  *x  •  V^N  k  •  V^N  =  S 
Dr  df  ^  * 


where  the  time  derivatives  x  and  k  are  related  to  the  frequency  u  through  the  characteristic 


equations 


i  =  V  . 


k  = 


Using  Cartesian  spacial  coordinates  (x.y),  and  polar  wavevector  coordinates  ik,d): 


(10a, b) 


0  =  tan'*  -2  , 


(11a, b) 


the  action  density  equation  (Eq.  9)  becomes 

aJV  ^  dN  dN  BN  BN  c 

where,  for  current  field  U  =  (m.v),  the  coefficients  c,-  are 


=  CgCOSO  +  u 


Cy  =  Cg  Sind  +  V , 


—  cos6  +  ^sina 


sinh(2kd)  I  Bx 


-  k  —cosH  +  (—  +  ^)cosasma  +  —  sin^al  , 
Bx  Bx  By  By  J 


sinh(2itrf)  ax 


Bd  .  a  Bd  n 
—  sina  -  —cos6 


^cos^O  +  (—  -  4^)cos^sin^  -  4^sin^^ 
By  By  Bx  Bx 
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where  we  have  used  Eqs.  1,  3,  4,  and  10  to  write  the  time  derivatives  of  x  and  k  in  terms  of 
the  group  velocity  and  gradients  of  the  depth  and  currents. 

2.5.  Source  Term 

The  source  term  S  includes  input  from  the  wind,  nonlinear  wave- wave  interactions, 
dissipation  from  such  processes  as  wave  breaking  and  bottom  friction,  surfactant  damping,  and 
other  processes  which  generate  or  dissipate  energy.  Although  the  different  processes  are 
inherently  interconnected,  following  Hasselmann  [1968],  the  source  is  often  written  as  a  sum  of 
separate  terms  (weak  interaction  approximation): 

S  =.  -  5^  *  5^  (17) 

where  5^  represents  input  from  the  wind,  describes  the  nonlinear  wave-wave  interactions, 
and  ^ds  includes  all  dissipation  terms.  Many  of  the  source  terms  are  poorly  understood  and  a 
number  of  competing  expressions  have  been  developed  to  describe  them.  The  terms  and  Sds 
are  usually  expressed  as  relatively  simple  analytic  functions  which  are  linear  in  N.  The 
nonlinear  term  ^nl  is  far  more  complicated  involving  integration  over  resonant  k-vectors  in  a  five 
dimensional  phase  space.  Perturbative  expressions  for  the  nonlinear  interaction  have  been 
derived  for  both  gravity  waves  [Hasselmann,  1962],  and  capillary  waves  [Valenzuela,  1972;  van 
Gastel,  1987].  However,  computing  these  integrals  directly  is  exceedingly  costly  and  becomes 
impractical  for  most  realistic  situations.  It  is  likely  that  little  progress  in  wind  wave  modeling 
will  be  made  with  these  first-principles  approaches  unless  significant  simplifications  to  the 
nonlinear  source  terms  can  be  found.  Presently,  a  few  such  simplifications  have  been  developed 
and  incorporation  of  these  schemes  into  our  model  is  in  progress. 

A  far  more  simple  form  for  the  source  term  which  is  often  sufficiently  accurate  for  the 
description  of  short  waves  in  mild  current  gradients  was  suggested  by  Hughes  [Hughes,  1978]. 
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In  this  ai>proach  the  source  is  expressed  as  a  Taylor  series  expansion  in  pow^  of  N.  To  second 
order  we  have 


5  *  a  +  /3iV  +  yN^. 

The  constant  term  a  has  been  associated  with  input  from  the  wind  through  the  Phillips  turbulent 
resonant  mechanism  (Phillips,  1957),  but  is  not  included  in  our  model.  Rather,  we  assume  that 
the  source  disappears  for  zero  action  density.  The  linear  term  fiN  is  then  the  dominant  term  Tor 
small  N  and  can  be  interpreted  as  the  initial  growth  term  due  to  input  from  the  wind.  The 
initial  gnnvth  of  waves  by  wind  has  been  investigated  in  detail  by  theory  and  experiments  (fo’ 
example  see  Plant,  [1982]).  Finally,  the  coefficient  y  can  be  set  by  the  requirement  that  the  total 
source  5  should  vanish  as  N  approaches  its  equilibrium  value  N^.  This  condition  fixes  7  to  - 
0/Ng.  The  Hughes  source  can  then  be  written: 


S  = 


-fiN 


(N  -  NJ 
No 


(19) 


Physically  this  source  tries  to  mimic  a  true  source  by  pushing  the  system  toward  a  known 
equilibrium  in  regions  of  uniform  current.  This  form  of  the  source  has  been  fully  implemented 
into  our  model.  However,  it  is  applicable  only  for  small  current  gradients  and  must  be  used 
with  caution  since  it  has  no  firm  foundation  in  first  principles. 

3.  Model 

To  solve  the  action  density  equation  we  use  a  finite  differencing  scheme  which  updates 
the  action  density  at  each  time  step.  The  action  density  equation  (12)  must  be  solved  for  a 
specified  grid  of  x,  y,  k,  6,  and  t.  Finite  step  methods  often  lead  to  instabilities  for  such 
multidimensional  systems  if  sufficient  care  is  not  exercised  in  the  numerical  procedure.  Our 
method  closely  follows  the  scheme  implemented  in  WAVEWATCH  described  by  Tolman 
[1991].  HowevCT,  we  use  the  {k,B)  space  formalism  rather  than  the  (w,0)  space  method 
implemented  in  WAVEWATCH.  In  this  scheme  the  propagation  terms  and  source  terms  are 
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handled  in  separate  time  steps.  The  propagation  is  done  using  a  predictor-corrector  Crank- 
Nickolson  procedure.  The  corrector  is  central  in  both  space  and  time  and  the  predictor  is 
forward  in  time  and  central  in  space.  To  illustrate  the  method  consider  the  simplified 
propagation  equation  (x-space  only): 

^  dN(x,t)  _  Q  (20) 

dt  *  dx 


The  action  is  updated  one  time  step  as  follows: 

Predictor: 


Half-Time  Value: 

,t„)  * 


(21) 


(22) 


Corrector: 


Arc 


(23) 


If  this  procedure  is  used  without  modification  spurious  solutions  or  negative  action  can 
occur.  As  discussed  in  Tolman  these  are  avoided  by  adding  a  diffusion  term  in  the  spacial 
dimensions,  using  a  limited  flux  approach  in  i^-space,  and  using  a  conservative  elimination 
scheme  in  0-space.  Adopting  the  more  compact  notation  N{ij,k,l,n)  =  V(x,  ,  yj  .  .  Bi ,  t„  ), 

the  complete  five-dimensional  action  density  equation  for  the  predictor  (the  corrector  has  the 
same  structure)  then  has  the  form: 
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where  a  is  the  upstream  fraction  (0  <  a  ^  1  for  c,  >  0  and  -1  <  a  <  0  for  c,  <  0),  and 
u  is  a  suffix  indicating  the  "upstream"  bin  in  ''-space  (next  bin  if  c*  >  0,  previous  bin  if  c*  < 
0),  The  sign  of  M  in  Eq.  30  is  the  sign  of  M*. 

The  first  order  upstream  scheme  used  in  the  propagation  in  x~space  assures  stable 
behavior  of  the  propagation  module  for  a  suitably  chosen  upstream  fraction  a.  The  stability  of 
the  propagation  is  governed  by  the  spacial  resolution  of  the  action  density  and  by  the  magnitude 
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of  the  source  term.  Since  the  diffusion  terms  reduce  the  accuracy  of  the  propagation  module 
it  is  best  to  keep  the  upstream  fraction  a  as  small  as  possible.  We  have  implemented  the 
expression  suggested  by  Tolman  [1991]  for  the  upstream  fraction: 


|a|  »  max[a^  ,  , 


=  min  1  ,  0. 10  +  3.75 


SJk,e)  Ax 
Fik,e)  c. 


where  is  chosen  large  enough  to  stabilize  the  source  free  propagation. 


The  propagation  in  k-space  incorporates  an  upstream  limited  flux  scheme.  In  this  method 
the  possibility  of  generating  negative  action  is  removed  by  requiring  that  the  flux  out  of  a  bin 
in  a  single  time  step  does  not  exceed  the  total  action  within  the  bin.  For  nearly  all  spectral  and 
spatial  bins  the  flux  limitation  will  normally  not  be  reached  and  M  and  M’  will  be  equal.  In  this 
case,  Eq.  (27)  reduces  to  the  normal  formulation  of  the  Crank-Nickolson  scheme.  Since  k  varies 
over  several  decades  we  have  chosen  an  exponential  k-  grid  for  better  sampling  of  the  lower 
wavenumbers  (ficquencies): 


k  =  k 

*  ‘  r 


{^) 


2t  \ 

\nax 


i^) 


it  =  1,2,  3,  ..,,N 


The  logarithm  of  this  ^-grid  is  uniformly  spaced.  Thus,  the  central  two  point  numerical 
approximation  is  appropriate  for  the  it-derivative  provided  we  write  it  as 

A  -i_?_  .  (33) 

dk  k  dln^ 

This  logarithmic  form  is  used  in  the  expression  for  (Eq.  27). 


The  propagating  in  0-space  (Eq.  28)  is  unaltered  from  the  normal  Crank-Nickolson 
derivative  formula.  To  stabilize  propagation  in  this  dimension  we  have  used  the  conservative 
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elimination  procedure.  This  amounts  to  setting  all  bins  with  negative  action  to  zero  while 
preserving  tlie  sum: 


2t 

N(Jc)  =  |iV(it,0)de.  (34) 

In  this  scheme  all  negative  action  for  a  given  J^-point  is  removed  first;  then  the  action  density 
for  all  directions  6  is  multiplied  by  a  constant  term  to  conserve  N{k). 

For  boundary  land  points  we  implemented  the  angle  derivative  upstream  scheme 
described  by  Tolman  [1991].  In  most  instances,  however,  we  have  been  working  in  patches  of 
the  ocean  where  no  land  points  exist.  In  these  cases  an  absorbing  boundary  is  used  where  the 
action  at  boundary  points  is  set  equal  to  the  action  at  adjacent  interior  points  at  each  time  step. 


At  the  boundaries  in  the  wavenumber  domain  we  assume  that  the  asymptotic  behavior 
of  the  spectrum  follows  that  of  an  assumed  equilibrium  spectrum  N^.  The  action  density  at  the 
upper  boundary  /r-points  is  generated  from  the  action  density  at  adjacent  interior  frequencies  by 

N{k)  =  NAk) ,  k>k^  (35) 

where  k^  is  a  user  supplied  upper  cutoff  wavenumber  which  depends  on  the  time  step  and 
resolution.  At  the  lowest  I:-point  where  the  action  density  is  rapidly  approaching  zero  the  action 
density  spectrum  is  simply  set  identically  to  the  equilibrium  spectrum. 

After  the  propagation  step  is  complete  the  source  term  is  included  by  integrating  the 
following  reduced  version  of  the  action  density  equation: 


10 


^  _ 

dt 


(36) 


The  integration  is  performed  using  a  simple  Euler  method  which  gives  the  updated  action  density 
^new  by  the  formula: 

(37) 

Here  is  the  action  density  output  from  the  propagation  module  (i.e.  without  the  source). 
Any  negative  action  that  results  from  the  integration  of  the  source  term  at  every  time  step  is 
simply  eliminated  by  setting  the  action  density  in  that  bin  to  zero.  The  time  step  is  chosen  small 
enough  to  satisfy  the  Courant  conditions,  At  <  Ax/Cf,  for  each  dimension,  {j:,}  =  {x,  y,  k,  0}. 
In  practice,  the  spectral  grid  is  chosen  fine  enough  so  that  the  upper  limit  on  At  is  constrained 
by  the  Courant  condition  in  the  spacial  domain. 

Updating  the  action  density  at  each  time  step  involves  a  call  to  the  propagation  module 
and  then  to  the  source  integrator  module.  The  procedure  is  illustrated  in  Fig.  1. 

4.  Testine 

The  Crank-Nickolson/predictor-corrector  propagation  scheme  has  been  thoroughly  tested 
previously  by  Tolman.  Our  code  was  ch«:ked  against  Tolman’s  WAVEWATCH  results  for 
each  of  the  compatible  test  cases  presented  in  his  thesis.  The  results  for  propagation  in  x-space 
and  fl-^ce  were  identical  in  appearance  to  Tolman’s  results.  The  propagation  in  it-space  could 
not  be  checked  against  WAVEWATCH  which  is  implemented  in  the  frequency  domain  rather 
than  wavenumber.  However,  internal  tests  of  the  I:-propagation  were  made  by  comparing  the 
limited  flux  scheme  to  several  other  forms  of  the  numeric  derivative. 

In  addition,  our  model  was  tested  against  the  time-independent  environmental  research 
institute  of  Michigan  (ERIM)  ocean  wave  model  (EOM).  The  EOM  solves  for  the  equilibrium 
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action  density  spectrum  for  a  given  steady  state  current  field.  The  comparisons  were  done  using 
the  one  dimensional  steady  current  field  of  the  form: 

u(x)  =  Ug  sech\x/S)  (38) 

which  is  a  stationary  soliton-like  current  which  goes  to  zero  at  the  boundaries  oo.  The 
parameters  and  8  were  chosen  to  be  0. 1  m/s  and  20  m.  A  plot  of  the  current  is  shown  in 
Fig.  2. 


To  estimate  the  changes  in  the  wind  waves  near  the  soliton  we  have  used  the  Hughes 
form  for  the  source  function  with  /3  given  by 


P  = 


0.04  2 

^AC 

fe-e  ] 

w 

K  U, 

vV/a 

[  2 

(39) 


where  u,  is  the  friction  velocity,  and  6^^  is  the  wind  direction  which  we  take  to  be  zero.  The 
equilibrium  spectrum  used  in  the  source  was  the  composite  spectrum  suggested  by  Bjerkaas  and 
Riedel  [1979]  which  is  given  in  appendix  A.  The  wind  direction  was  along  +x  with  a  speed  of 
5  m/s.  A  uniform  95  point  grid  with  Ax  =  4  m  was  used  for  the  mesh  along  x.  The  radial 
wavenumber  mesh  is  computed  with  Eq.  (32)  using  82  Jt-points,  with  X,„i„=  0.(X)1  m,  and 
\nax=50  m.  The  angular  spectral  mesh  is  uniform  with  A0=1O°  (36  angles). 


Before  the  action  density  equation  can  be  solved  it  is  necessary  to  provide  appropriate 
boundary  conditions  for  each  of  the  five  dimensions.  The  simulation  was  started  by  applying 
the  equilibrium  spectrum  to  all  95  points  along  x,  although  any  other  initial  conditions  should 
in  principle  produce  the  same  equilibrium  result.  On  subsequent  time  steps  only  the  boundary 
x-points  were  set  by  applying  an  absorbing  boundary  condition  as  discussed  previously.  The 
boundaries  in  the  spectral  domain  were  set  according  to  Eq.  (35)  with  =  2t/(0.005  m). 

To  examine  changes  in  the  action  density  we  have  defined  a  fractional  difference 
spectrum  R  as  follows: 
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(40) 


R(^,k,e) 
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where  is  the  wind  direction.  Ilie  time  development  of  /?  as  a  function  of  position  x  is  shown 
for  wavelengths  ranging  from  1  cm  to  15  m  in  Fig.  3.  The  solid  line  in  each  plot  is  the  function 
R  at  t  le  r=2S0  s,  while  the  dotted  lines  are  snapshots  of  R  at  previous  times  separated  by  50  s 
intervals.  Apparent  in  the  plots  are  several  transients  which  have  mostly  propagated  out  of  the 
plot  window  by  the  maximum  time  (250  s).  Equilibration  rates  are  a  function  of  wavelength, 
but  all  wavelengths  appear  to  have  reached  equilibrium  by  250  s.  We  find  that  wavelengths  near 
the  center  of  the  spectrum,  X  *  1  m  to  2  m,  are  most  strongly  perturbed  by  the  current  field. 
The  functional  form  of  the  plots  are  consistent  with  the  rule  that  longer  waves  are  governed  by 
the  current  amplitude,  while  shorter  waves  are  governed  by  the  current  gradient.  These  limiting 
forms  are  easily  derived  from  the  action  density  equation  (9,12,  also  see  Lyzenga,  [1991]), 
which  for  small  currents,  and  current  gradients  (and  thus  for  /?<<!)  leads  to  the 
approximations 

►  Cjj  ►  (short  wave  refractive  limit) , 

"  (41) 

/35 ,  u  <  c„  Gong  wave  advective  limit) , 

o* 


R  - 
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where  ii  is  the  dimensionless  number, 


v(x,k,B) 
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(42) 


which  is  negative  for  wavelengths  anywhere  in  the  tail  of  the  spectrum. 


A  comparison  between  the  NTOW  model  results  and  those  of  the  EOM  are  shown  in 
Fig.  4.  The  NTOW  model  results  are  plotted  at  time  t=250  s  when  the  system  appears  to  have 
reached  equilibrium  near  x=0.  The  comparison  between  the  models  is  quite  good  considering 
the  differences  in  the  numerical  methods.  The  positive  "bump"  which  appears  in  the  NTOW 
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results  for  1  m  waves  and  at  the  edge  of  the  plot  for  2  m  waves  are  transients  which  are 
propagating  out  of  the  plots.  The  small  differences  between  the  results  of  the  two  models  for 
the  shortest  wavelengths  are  due  to  differences  in  the  way  the  tail  region  of  the  ^)ectrum  is 
handled. 

Figures  Sa,  b,  and  c  show  the  full  spectrum  results  at  the  center  of  the  soliton  and  at 
X  —  T  20  m  for  (5a)  R  using  the  NTOW  model,  (5b)  R  using  the  EOM,  and  (5c)  the  difference 
between  R  for  the  two  models.  It  is  clear  from  these  plots  that  both  models  give  very  similar 
results  for  the  fractional  difference  spectrum.  However,  the  maximum  perturbation  of  the 
spectrum  for  EOM  is  slightly  larger  than  that  for  the  NTOW  model.  Still,  from  Fig.  5c  it 
appears  that  the  maximum  fractional  difference  between  the  two  models  is  less  than  5  %  for  all 
regions  of  the  spectrum. 

A  radar  image  of  the  soliton  feature  was  derived  from  the  wave  spectrum  using  the 
composite  scattering  model  of  Wright  [1968],  which  is  implemented  in  a  module  of  the  EOM. 
Radar  imagery  can  be  derived  using  the  output  wave  field  from  either  wave  propagation  model 
as  input  to  the  EOM  radar  module.  Figure  6  compares  the  RAR  cross  section  per  unit  area 
obtained  with  the  two  ocean  wave  models  for  X-Band,  HH-polarization  and  an  incidence  angle 
of  37“.  We  find  that  the  maximum  cross  section  variation  for  the  EOM  is  larger  than  that  for 
the  NTOW  model  consistent  with  the  larger  perturbations  in  the  wave  action  density  found  for 
EOM.  Nevertheless,  the  two  models  give  very  similar  results  for  the  radar  return,  especially 
in  the  general  shape  of  the  curve. 

The  results  for  the  stationary  soliton  show  good  agreement  between  the  EOM  and  the 
NTOW  model.  This  suggests  that  the  NTOW  model  should  give  satisfactory  results  in  cases 
where  the  currents  or  winds  are  not  steady.  In  the  next  section  we  examine  such  a  case  for  a 
current  "Rip"  formed  by  the  subduction  of  Gulf  Stream  water  under  shelf  water. 
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S.  Rip  Feature 


5.1.  Background 

During  the  NRL  high-resolution  remote  sensing  experiment  real  aperture  radar  (RAR) 
data  was  collected  over  the  shelf  region  off  Cape  Hatteras.  The  RAR  data  show  many  features 
which  are  known  to  be  associated  with  current  structures  near  the  Gulf  Stream.  However,  the 
precise  nature  of  these  features  is  still  unknown.  Some  of  the  more  obvious  mechanisms  which 
could  produce  these  features  are,  wave  current  interaction,  wave  breaking,  and  surfactant 
damping  or  combinations  of  these  effects. 

A  particularly  striking  feature  seen  in  the  NRL  high-resolution  RAR  data  was  a 
meandering  line  in  the  region  between  the  Gulf  Stream  water  and  the  coastal  shelf  water.  The 
RAR  image  of  the  line  feature,  which  has  been  named  a  current  "Rip"  feature,  is  shown  in 
Fig.  7.  The  Rip  is  from  5dB  to  20dB  above  the  background  and  is  visible  over  several 
kilometers  along  the  shelf.  Unfortunately,  no  accurate  measurements  of  surface  currents  are 
available  in  this  region.  From  measurements  of  deeper  ocean  currents  and  analysis  of  the  RAR 
data  it  was  estimated  that  the  Rip  is  advected  primarily  in  a  northerly  direction  at  about  0.4  m/s. 
It  is  believed  that  the  Rip  feature  results  from  wave-current  interaction  effects  associated  with 
surface  convergence  (downwelling)  «'org  a  line  where  Gulf  Stream  water  is  subducted  under 
slower  shelf  water  [Mied  et  al.,  1992].  The  subduciion  occurs  because  the  Gulf  Stream  water 
is  more  dense  than  the  shelf  water  due  to  its  higher  salinity. 

5.2.  Modeling  of  Rip 

Although  the  RAR  intensity  has  some  variation  along  the  Rip,  the  effect  is  expected  to 
be  predominantly  associated  with  current  variat!-.*  ^  i"  *he  direction  perpendicular  to  the  feature. 
We  have,  therefore,  modeled  the  Rip  currents  as  one  dimensional.  The  form  used  to  describe 
the  currents  in  the  direction  perpendicular  to  the  Rip  is 
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v(y)  ~  a  -  b  tanh 


y-ct 

5 


(43) 


A  plot  of  the  current  is  shown  in  Fig.  8.  The  parameters  a,  b.  c,  and  5  were  estimated  to  be 
near  the  values  a =0.5  m/s,  h  =  0.1  m/s,  c  =  0.4  m/s,  and  6  =  10.  m.  However,  these  values 
are  only  approximate  since  the  available  current  data  near  the  rip  is  inadequate  to  determine  the 
current  field  precisely.  For  this  reason,  we  have  examined  the  sensitivity  of  the  radar  return  to 
variations  in  these  parameters. 

The  calculations  were  done  using  the  NTOW  model  in  a  coordinate  system  with  x 
pointing  along  the  Rip  (predominantly  ,  east)  and  y  perpendicular  to  the  Rip  (predominantly 
north).  The  grid  in  the  spacial  domain  along  y  uses  145  points  with  a  uniform  spacing  of 
Ay  =  4  m.  The  radial  wavenumber  mesh  is  computed  from  Eq.  (32)  using  62  it-points,  with 
Xn,in=  0.001  m,  and  X^=  50  m.  The  angular  spectral  mesh  is  uniform  with  24  angles 
=  15®).  The  simulation  was  initiated  by  applying  the  equilibrium  spectrum  to  all  points 
at  time  /  =  0.  On  subsequent  time  steps  was  applied  only  to  boundary  points. 

To  estimate  changes  in  the  wind  waves  near  the  Rip  we  have  again  used  the  Hughes  form 
with  /3  given  by  Eq.  (39)  and  the  Bjerkaas  and  Riedel  [1979]  equilibrium  spectrum  (see 
appendix  A).  The  measured  wind  direction  and  speed  are  shown  in  Fig.  9.  From  Fig.  9  it 
follows  that  over  typical  equilibration  time  scales  for  the  Bragg  waves  and  waves  of  a  few 
meters  the  wind  variation  is  negligible.  At  the  time  the  RAR  data  were  obtained  the  wind  was 
from  the  South-West  at  about  3  m/s.  In  the  calculations  we  use  a  fixed  wind  speed  of  3  m/s  and 
have  also  ignored  the  relative  wind  speed  corrections  in  the  source  term  since  the  current  speed 
is  small. 

Some  preliminary  calculations  for  the  Rip  were  also  done  using  the  EOM  wave 
propagation  model  by  shifting  the  coordinate  system  to  the  frame  centered  on  the  rip.  In  this 
frame  the  currents  are  steady  so  that  the  EOM  model  is  applicable.  However,  our  calculations 
using  the  EOM  show  very  poor  convergence  of  the  perturbed  wave  field  with  grid  spacing  for 
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many  of  the  cases  considered  here.  Although  the  shape  of  the  wave  field  predicted  with  the 
EOM  is  consistent  with  NTOW  results,  the  magnitude  continued  to  change  significantly  for 
variations  in  grid  spacing,  even  when  the  grid  spacing  was  much  less  than  the  convergence  zone 
width.  Due  to  its  slow  convergence,  the  wave  field  predictions  from  the  EOM  model  are  not 
included  in  this  report.  Further  discussion  of  this  problem  is  given  in  appendix  B. 

5.3.  Wave  Field  Prediction 

Figure  10  shows  the  time  development  of  the  y-component  of  the  fractional  difference 
spectrum,  R(k,d=9QP),  for  a  number  of  wavelengths  ranging  from  1  cm  to  10  m.  The  wind  is 
at  45°  toward  the  north  east  (i.e.  6^=45°).  The  different  curves  for  a  given  wavelength 
represent  the  function  R  at  times  separated  by  A/  =  30  s.  The  shift  in  the  curves  toward  the 
north  (higher  y)  is  due  to  the  advection  of  the  Rip  with  a  drift  speed  of  c  =  0.4  m/s.  The 
general  behavior  of  R  as  a  function  of  wavelength  can  be  understood  from  the  approximate 
expressions  given  in  Eq.  (42).  That  is,  we  find  that  the  shortest  waves  are  governed  by  the 
magnitude  of  the  current  gradient  while  the  longest  waves  are  governed  by  the  current 
magnitude.  As  the  waves  approach  the  Rip  from  the  south  they  are  refracted  towards  the  north 
by  the  current  gradient  according  to  Eq.  (10).  Longer  waves  are  also  shifted  to  higher 
frequencies  upon  entering  the  slower  moving  water  on  the  north  side.  Thus,  shorter  waves  are 
enhanced  near  the  Rip  edge  on  the  north  side.  The  results  also  show  that  the  shortest  waves  are 
affected  only  very  near  the  vicinity  of  the  convergence  front  by  the  current  gradient,  while 
longer  waves  are  advected  and  continue  to  propagate  with  a  perturbed  action  density  some 
distance  from  the  Rip  on  the  north  side. 

Even  for  a  fixed  wind  direction  the  angle  that  the  wind  makes  with  the  Rip  varies  since 
the  Rip  meanders.  We  have  therefore  computed  the  action  density  across  the  Rip  for  several 
wind  directions.  Figure  1 1  shows  the  contour  and  surface  plots  of  R  for  wind  directions  towards 
(1  la)  20°,  (11b)  45°,  and  (1  Ic)  70°  from  the  Rip  line  (x-axis).  For  each  wind  direction  the  plots 
are  shown  at  locations  y  =  -20  m  (south  side  of  Rip),  y  =  0  m  (center  of  Rip),  and  y  =  +20 
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m  (north  side  of  Rip).  The  largest  enhancement  of  the  waves  is  on  the  north  side  and  is  due  to 
the  shift  of  the  waves  to  shorter  wavelengths  upon  crossing  the  convergence  zone  center.  This 
oihancement  will  be  greatest  for  waves  traveling  perpendicular  to  the  Rip  line.  Thus,  the  peak 
on  the  north  side  is  found  to  be  largest  for  the  wind  direction  6^  =  70°. 

The  sensitivity  of  R  to  variations  in  the  convergence  speed  b  or  width  5  is  presented  in 
Fig.  12.  The  spectrum  R  is  shown  for  the  matrix  of  parameters  b  =  .1  m/s,  .2  m/s,  and  .3  m/s, 
and  6  =  10  m,  20  m,  and  30  m  for  fixed  a  (=.5  m/s),  c  (=.4  m/s),  and  (=  45°).  In  all 
cases  the  general  shape  of  the  fractional  difference  spectrum  is  very  similar  to  Fig.  11b. 
However,  the  magnitude  of  the  enhancement  varies  with  these  parameter  changes,  increasing  as 
the  convergence  speed,  b,  increases  and  decreasing  as  the  convergence  zone  width  5  increases. 
The  sensitivity  of  R  to  variations  in  b  appears  to  be  especially  strong  and  nonlinear,  growing  by 
threefold,  for  example,  for  a  doubling  of  b  from  .  1  m/s  to  .2  m/s.  However,  it  is  expected  that 
for  such  large  current  gradients  the  Hughes  source  term  may  be  inappropriate.  The  results  for 
very  large  fractional  changes  should  therefore  be  viewed  with  some  caution. 

5.4.  Radar  Predictions 

We  have  also  derived  composite  scattering  model  predictions  of  the  RAR  imagery  for 
each  of  the  cases  involving  the  Rip  discussed  above.  In  these  calculations  the  plane  flies  parallel 
to  the  Rip  on  the  north  side  as  was  done  in  the  measured  RAR  data  shown  in  Fig.  7.  The  RAR 
incidence  angle  at  the  Rip  was  near  37°,  but  we  have  examined  the  response  of  the  RAR  to 
variations  in  this  angle.  The  predicted  HH-polarized  X-band  (12.7  GHz)  RAR  return  for  the 
Rip  when  the  standard  Rip  parameters  (a=.5  m/s,  b=.\  m/s,  c=.4  m/s,  6  =  10  m,  =  45°) 
and  a  37°  incidence  angle  are  assumed  is  shown  in  Fig.  13.  Consistent  with  the  simplest  radar 
models,  we  find  that  the  shape  of  the  curve  looks  much  like  the  shape  of  R  for  short  waves. 

Figure  14  shows  the  dependence  of  the  RAR  on  wind  direction  relative  to  the  Rip  (14a), 
convergence  speed  b  for  fixed  S  =  10  m  (14b),  and  RAR  incidence  angle  (14c).  It  is  found  that 
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the  changes  in  the  RAR  cross  section  roughly  follow  changes  in  R  as  the  wind  direction  and 
convergence  speed  are  varied  (14a, b).  As  in  the  behavior  of  R,  the  RAR  cross  section  is 
{Mirticularly  sensitive  to  changes  in  the  convergence  speed,  growing  nonlinearly  with  increases 
in  parameter  b  (Fig.  14b).  The  largest  return  is  for  convergence  speed  b=.3  m/s  (with  6  =  10 
m)  where  the  gain  is  about  12  dB.  This  is  well  within  the  range  of  the  measured  gain  of  S  dB 
to  20  dB.  The  dependence  on  incidence  angle  shows  a  downward  trend  in  the  return  near  37° 
incidence.  The  falloff  in  the  curve  in  Fig.  13c  near  zero  incidence  is  due  to  a  breakdown  in  the 
composite  scattering  model  for  very  small  incidence  angles  and  is  not  realistic. 

A  more  complete  view  of  the  dependence  of  the  HH-RAR  return  on  variations  of 
parameters  b,  and  6  is  shown  in  Fig.  IS.  Here  we  show  the  maximum  gain  in  the  RAR  return 
for  the  3X3  matrix  of  cases  with  different  b  and  5  values  and  a  fixed  incidence  angle  of  37"  and 
wind  direction  of  45°.  As  expected  the  largest  gain  in  cross  section  (about  12  dB)  is  obtained 
for  the  smallest  5  and  largest  b  considered.  The  dependence  on  b  is  strong  as  discussed 
previously.  However,  the  sensitivity  of  the  RAR  to  changes  in  6  is  relatively  weak.  In  general, 
our  results  for  the  maximum  radar  return  in  the  Rip  appear  to  be  somewhat  less  than  those 
observed  in  the  RAR  data.  This  is  not  surprising  since,  among  other  factors,  we  have  neglected 
wave  breaking,  which  is  expected  to  be  important  in  the  Rip  region,  especially  for  the  higher 
convergence  speeds.  Surfactant  effects  are  also  not  included  in  our  simulations,  however, 
preliminary  studies  of  expected  surfactant  effects  suggest  these  will  be  weak  near  the  Rip,  but 
could  be  significant  away  from  the  Rip,  reducing  the  amplitude  of  waves  shorter  than  about  1  m. 

6.  Summary 

A  time  dependent  ocean  wave  (NTOW)  model  has  been  developed  to  simulate  the 
propagation  of  waves  over  current  and  depth  structures.  The  model  solves  the  action  density 
conservation  equation  using  a  finite  step  method.  Comparisons  of  the  model  were  made  against 
Tolman’s  WAVEWATCH  model  and  the  ERIM  Ocean  Model  (EOM)  to  determine  the  reliability 
of  the  propagation  scheme.  Agreement  between  the  various  models  was  found  to  be  satisfactory. 
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We  have  used  the  NTOW  model  along  with  the  composite  scattering  model  to  make  predictions 
of  the  radar  return  from  a  current  Rip  and  compared  our  results  with  measured  RAR  data.  Our 
predictions  of  the  RAR  cross  section  for  the  current  Rip  are  found  to  be  consistent  with 
measurement,  but  generally  show  smaller  variation.  However,  the  effects  of  whitecapping  and 
surfactants  have  not  been  included  in  our  modeling.  It  is  expected  that  including  these  effects 
will  produce  better  agreement  with  experiment.  It  is  hoped  that  these  effects  can  be  included 
in  ^ture  work.  Also,  more  realistic  source  terms  are  being  implemented. 
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APPENDIX  A:  Equilibrium  Spectrum 


The  equilibrium  action  spectral  density  Nf,(k,6)  is  related  to  the  equilibrium  variance 
density  spectrum  FJik,6)  by 

NJ,k,V)  .  (AD 

The  spectrum  FJk.d)  is  that  proposed  by  Bjerkaas  and  Riedel  [1979]  which  for  low  or  moderate 
wind  speeds  <  10  m/s)  is  written  as: 

=£.(*)  (A2) 


where 


G,(d) 


and 


Eo(k) 


>  I:  < 

E^k),  k,<k<k, 

£v(A:) ,  k^k^ 


where  Epj^k)  is  the  Pierson-Moskowitz  spectrum: 


Ep^k)  = 


0.0081 
2k* 


-exp 


0.74g^ 
{30  u,)*k^ 


(A3) 


(A4) 


(A5) 
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E/^k)  is  the  Kitaigorodskii  spectrum: 


where  k„  is  taken  to  be  i^=363  m  ‘,  and  p=S-Logxdu*. 
Eyifii)  is  the  Mitsuyasu-Honda  spectrum: 


and  Ey(k)  is  the  viscous  damping  spectrum 


£v(k)  = 


0.0001473  («,  klf 


The  values  of  ki  ,  k^  ,  and  k^  satisfy  the  matching  conditions: 


(A6) 


(AT) 


(A8) 


EmM  =  ^v(^)-  (A9) 

This  fixes  ki  and  ^3  to  values  which  are  functions  of  u».  However,  ^2  is  arbitrary  and  was  set 
to  the  value  ^2=263.86  m.,.  The  friction  velocity  u,  is  related  to  the  wind  speed  at  height 
h  above  the  sea  surface  through  the  formula 

(AlO) 


where 


(All) 


-  =  6.84  X  IQ-^  ^  ^  28  X  10-3  ^2  _  ^  ^ 

u. 

Typically,  h  is  chosen  to  be  19.5  m  in  which  case  this  reduces  approximately  to  (for  ~  S 
m/s) 

«^(/i=19.5m)  *  30  u,  .  (All) 

A  plot  of  the  equilibrium  variance  density  spectrum  Foik,d)  is  shown  in  Fig.  Al. 
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APPENDIX  B:  Convergence  of  Wave  Models  with  Grid  Spacing 


Wave  models  of  a  convergent  flow  with  current  ^)eeds  near  that  of  the  Rip  frequently 
predict  very  large  enhancements  in  the  wave  field  when  a  Plant-Hughes  source  is  used.  We 
have  used  two  models  to  make  these  predictions:  the  ERIM  Ocean-wave  Model  (EOM) 
RSLWAVE,  and  the  NRL  Time  dependent  Ocean  Wave  (NTOW)  model.  Our  calculations  have 
consistently  found  some  disagreement  between  these  two  models,  especially  for  the  larger 
convergence  speeds.  RSLWAVE  always  produces  larger  enhancements  than  NTOW,  but  the 
differences  between  the  models  are  negligible  when  the  fractional  difference  spectrum,  R,  is 
small: 


R(x,k) 


N(x,k)  -  N^(k) 

N,ik) 


<<  1, 


(Bl) 


where  is  the  equilibrium  spectrum.  The  disagreement  between  the  models  grows 
substantially,  however,  as  R  increases  (for  Ax  =  4  m).  For  example,  for  the  largest  Rip 
velocity  gradients  considered,  0.03  s‘‘,  RSLWAVE  predicts  an  R  value  of  around  16,  while 
NTOW  predicts  10.2  -  a  dramatic  difference! 


From  some  preliminary  testing,  the  source  of  this  discrepancy  between  the  models  is 
apparently  due  to  the  high  sensitivity  of  RSLWAVE  to  the  choice  of  grid  spacing  Ax.  For  any 
realistic  modeling  the  grid  spacing  should  be  kept  small  compared  to  the  width  of  the 
convergence  zone  5.  But,  it  is  found  that  the  predictions  for  R  using  RSLWAVE  vary 
appreciably  with  changes  in  Ax  even  when  this  condition  is  well  satisfied.  Figure  Bl  shows  a 
plot  of  the  peak  value  of  R  as  a  function  of  5/ Ax  at  a  point  10  m  on  the  north  side  of  the  Rip. 
The  Rip  current  is  modeled  as 

«(x)  =  0.1  -  0.3  tanh(x/5) 

with  6  =  10  m.  We  see  that  convergence  in  the  NTOW  model  is  reached  around  Ax  =  4  m, 
whereas  RSLWAVE  is  still  changing  significantly  in  this  region. 


29 


The  reason  for  the  slow  convergence  of  RSLWAVE  is  not  related  to  round  off  errors. 
Doubling  the  precision  made  no  change  in  the  results.  The  cause  must  therefore  be  related  to 
the  propagation  scheme  used  in  the  models.  RSLWAVE  uses  a  forward  propagation  scheme  in 
which  spectral  derivatives  of  the  action  density  are  connected  in  the  discrete  relation  to  spacial 
derivatives  at  different  points  in  phase  space.  Thus  a  correlation  between  spacial  resolution  and 
spectra]  variations  is  made  possible,  and  should  grow  worse  with  increasing  Ax.  In  NTOW  all 
spectral  and  spacial  derivatives  are  central  derivatives  about  the  same  point  in  phase  space  in  the 
discrete  expression  of  *he  action  density  equation.  This  minimizes  the  possibility  of  any 
correlation  between  the  spectral  changes  and  the  grid  resolutions.  In  general  RSLWAVE  tends 
to  produce  larger  perturbations  in  the  wave  field  due  to  current  structures  even  when  well 
convayed.  This  may  be  related  to  die  use  of  a  diffusion  term  in  NTOW  which  tends  to  smooth 
rapid  variations  in  the  action. 

Of  course,  for  the  cases  considered  here  it  should  be  realized  that  the  Plant-Hughes  form 
may  be  inappropriate.  In  any  case,  a  user  of  these  models  should  be  aware  of  these 
discrepancies. 
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Fig.  1.  Time  loop  diagram  for  the  NTOW  model. 


Fig.  3.  Time  development  of  fractional  difference  spectrum  vs  x  for  several 

wavelengths  (WL  values).  The  solid  line  is  R  at  time  r-2S0  s.  Dotted  lines  are  plots  of  R  for 
previous  times  separated  by  50  s. 
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Fig.  4.  Comparison  between  EOM  (dotted  lines)  and  the  NTOW  model  (solid  lines)  results  for 
the  fractional  difference  spectrum  vs  x  at  several  wavelengtfis  (WL’s). 
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Fig.  5a.  NTOW  model  contour  and  surface  plots  of  equilibrium  fractional  difference  spectrum 
Rikjd)  for  the  soliton  at  a:  -  -20  m,  a:  -  0  m  (center  of  soliton),  and  at  -  20  m.  The  axes  are  the 
dimensicmless  numbers  Kx  -  ln(t/it„^os0,  and  Ky  -  ln(it/ifc^sin0,  where  k^  -  (2jt/.(X)l)  m‘‘. 
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x= 


-20  m 


X=  0  m 
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Fig.  5b.  EOM  contour  and  surface  plots  of  fractional  difference  spectrum  R(k,0)  for  the  soliton 
at  x  -  -20  m,  j:  ->  0  m  (center  of  soliton),  and  jr  -  20  m.  The  axes  are  the  dimensionless  numbers 
Kx  -  ln(i{://i;„^os0,  and  Ky  -  ln(ifc/iL..)sin6.  where  -  (2to'.001)  m  *. 
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Fig.  5c.  Contour  and  surface  plots  of  the  spectrum  (A^eom  •  ^ntowV^^*  at  equilibrium  for  the 
soliton  at  ;c  -  -20  m,  -  0  m  (center  of  solitcm),  and  x  -  20  m.  The  axes  are  the  dimensionless 
numbers  Kx  -  ln(fcJi„^cos0,  and  Ky  -  ln(ik/i^^sin0,  where  -  (2n/.(X)l)  m  ‘. 
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(S/UJ)  (fc)A 


Current 


y  (nr.) 


Fig.  8.  Initial  Rip  current  for  standard  parameters:  a  -  .5  m/s,  b  ~  A  m/s,  5  -  10  m  in  Eq.  43. 
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WIND  MEASUREMENTS  NEAR  RIP  FEATURE 
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Fig.  9.  Wind  speed  and  direction  near  Rip  on  day  RAR  data  was  collected. 
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Fig.  10.  Development  of  fractional  difference  spectrum  /?(^,9-90°)  for  1  cm,  50  cm,  1  m,  S  m, 
7  m,  and  10  m,  waves  propagating  across  a  current  "Rip"  (parameters  a  -  0.5  m/s,  £>  -  0.1  m/s, 
c  -  0.4  m/s,  and  8-10  m).  The  wind  direction  is  e„  -  45°  toward  ME.  The  time  interval 
between  the  curves  on  each  plot  is  A/  -  30  s.  Notice  diat  the  effect  m  the  shortest  waves  is 
small  and  very  localized  to  the  Rip  regitm.  For  medium  size  waves  (-Im)  the  effect  is  larger 
and  less  localized.  For  still  longer  waves  the  effect  is  reduced  but  extends  well  into  the  north 
side  region  (y  >  0)  of  the  Rip. 
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Fig.  11a.  NTOW  model  contour  and  surface  plots  of  equilibrium  fractional  difference  spectrum 
for  3  m/s  winds  blowing  toward  0^-20°  across  the  Rip.  R  is  shown  at  y  -  -20  m  (south 
side),  -  0  m  (center),  and  -  20  m  (north  side).  The  axes  are  the  dimensionless  numbers  Kx 
-  In(i/^„)cos0,  and  Ky  -  ln(^^„)sin0,  where  -  (27i/.001)  m  *.  The  Rip  parameters  are 
a  -  .5  ni/s,  A  m/s,  c  -  .4  m/s,  5  -  10  m. 
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Fig.  13.  X-Band  HH-RAR  cross  section  per  unit  area  at  3T  incidence  for  standard  Rip 
parameters:  a  -  .5  m/s,  b  ^  A  m/s,  c  -  .4  m/s,  5  -  10  m.  The  flight  was  parallel  to  the  Rip  on 
the  north  side.  The  wind  was  3  m/s  from  45*’  south  of  west 
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Fig.  14a.  Dependence  of  X-Band  HH-RAR  return  oi  wind  dinectitxi  for  Rip  feature. 
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RIP:  Dependence  on  Convergence  Veto 
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Fig.  14b.  Dependence  of  X-Band  HH-RAR  return  on  convergence  speed  for  Rip  feature. 
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Cross  Section  Gain  (dB) 


Fig.  14c.  Dependence  of  X-Band  HH-RAR  return  on  incidence  angle  for  Rip  feature. 
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Fig.  15.  Surface  plot  showing  the  maximum  X-Band  HH-RAR  cross  section  gain  (above 
background)  for  variatitms  in  the  ctmvergence  speed  (parameter  b),  and  convergence  zone  width 
(parameter  5).  All  other  Rip  parameters  were  set  to  the  standard  values. 
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Fig.  Al.  Equ_ 
speed  of  5  m/s. 
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Fig.  Bl.  Peak  value  in  the  fractional  difference  spectrum,  R^,  at  x  -  10  m,  predicted  by  the 
models  RSLWAVE  and  NTOW  as  a  function  of  grid  spacing  Ax.  The  Rip  current  model  is  used 
with  the  parameters  a  -  0.1  m/s,  b  -  0.3  m/s,  c  -  0.  ^s,  and  6  -  10  m.  The  wind  is  at  3  m/s 
from  45®  south  of  west  The  peak  corresponds  to  a  wavelength  of  about  1.2  m. 
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